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^ ■ Abstract: We study the estimation of a stable Cox-Ingersoll-Ross model, which 

is a special subcritical continuous-state branching process with immigration. The 
process is characterized in terms of some stochastic equations. The exponential 
ergodicity and strong mixing property of the process and the heavy tail behavior 
of some related random sequences are studied. We also establish the convergence 
of some point processes and partial sums associated with the model. From those 
fV| . results, we derive the consistency and central limit theorems of the conditional 

Oh ! least squares estimators and the weighted conditional least squares estimators of 

^ | the drift parameters based on low frequency observations. A weakly consistent 

estimator is also proposed for the volatility coefficient based on high frequency 
observations. 
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X ! 1 Introduction 



The Cox-Ingersoll-Ross model (CIR-model) introduced by Cox et al. (1985) has been used 
widely in the financial world. This model has many appealing advantages. In particular, 
it is mean-reverting and remains positive. Let a > 0, b > and a > be given constants. 
The classical CIR-model is a positive diffusion process {X(t) : t > 0} defined by 

dX(t) = (a- bX(t))dt + ay/X(t)dB(t), (1.1) 

where {B(t) : t > 0} is a standard Brownian motion. The process defined by (II. ip has 
continuous sample paths and light tailed marginal distributions. 
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It is well-known that many financial processes exhibit discontinuous sample paths and 
heavy tailed distributions. Let (Q, J^, J^, P) be a filtered probability space satisfying the 
usual hypotheses. A natural generalization of (11. ip is the stochastic differential equation 

dX t = (a- bX t )dt + a ^/X^ldZ u (1.2) 

where {Z t : t > 0} is a spectrally positive stable (J^)-Levy process with index 1 < a < 2. 
For a = 2, we understand the noise as a standard Brownian motion, so (II. 2p reduces to 
(11.11) . When 1 < a < 2, we assume it is a stable process with Levy measure 

"«W ■= a v { -a)z^ (L3) 

By a result of Fu and Li (2010), there is a pathwise unique positive strong solution 
{X t : t > 0} to (11.21) . We refer to this process as a stable Cox-Ingersoll-Ross model 
(SCIR-model). We shall see that the discontinuous SCIR- model indeed captures the 
important heavy tail property. The reader may refer to Borkovec and Kliippelberg (1998), 
Embrechts et al. (1997, Section 7.6) and Fasen et al. (2006) for similar modifications of 
the CIR-model. The SCIR-model is a particular form of the so-called continuous- state 
branching processes with immigration (CBI-processes), which arise as scaling limits of 
Galton-Watson branching processes with immigration (GWI processes); see, e.g., Kawazu 
and Watanabe (1971). The general CBI-processes were also constructed and studied in 
terms of stochastic integral equations in Dawson and Li (2006, 2012), Fu and Li (2010) 
and Li and Ma (2008). 

The estimation for stochastic processes based on the minimization of a sum of squared 
deviations about conditional expectations was developed in Klimko and Nelson (1978). 
They applied their results to the conditional least squares estimators (CLSEs) of the 
offspring and immigration means of subcritical GWI processes. Their estimators are es- 
sentially the same as those studied by Quine (1976, 1977). By the results of Klimko 
and Nelson (1978) and Quine (1976, 1977), under a finite third moment condition, as 
the sample size n goes to infinity, the errors of the CLSEs decay at rate n~ x l 2 and they 
are asymptotically Gaussian; see also the earlier work of Heyde and Seneta (1972, 1974). 
The asymptotic properties of CLSEs of GWI processes with general offspring laws were 
studied in Venkataraman (1982) and Wei and Winnicki (1989). Based on the idea of Nel- 
son (1980), the weighted conditional least squares estimators (WCLSEs) of the offspring 
and immigration means of GWI processes were proposed by Wei and Winnicki (1990), 
who proved some self-normalized central limit theorems for the estimators. The limiting 
distributions in Wei and Winnicki (1990) are also Gaussian except in the critical case. 
The reader can refer to de la Pena et al. (2009) for recent developments on self-normalized 
limit theorems and their statistical applications. The estimation problems of the CIR- 
model defined by (II. ip were studied by Overbeck and Ryden (1997). They proposed some 
CLSEs and WCLSEs and proved a Gaussian central limit theorem for them; see also 
Overbeck (1998). 

In this work, we give some estimation of the drift coefficients (b, a) of the SCIR-model 
using low frequency observations at equidistant time points {kA : k — 0, 1, • ■ • , n} of a 
single realization {X t : t > 0}. For simplicity, we take A = 1, but all the results presented 
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below can be modified to the general case. We shall also consider the parameters 

l = e~\ p = ab~ 1 (l- 1 ). (1.4) 

Following Klimko and Nelson (1978) and Overbeck and Ryden (1997), we first define the 
CLSEs of the parameters. The basic ideas are explained as follows. By applying Ito's 
formula to (11. 2p . for any t > r > we have 



(1.5) 



X t = e-W-^Xr + a f e'^ds + a f e'^X^dZ^ 

J r J r 

From (11.51) we obtain the stochastic regressive equation 

X k = p + 7X fe _i + e k , (1.6) 

where 

e k = a I" e- b ^XlL a dZ s . (1.7) 



One can see that {e k '■ k > 0} is a sequence of martingale differences with respect to 
{J^fc : k > 0}. The CLSEs of (7, p) and (b, a) can be given by minimizing the sum of 
squares 



J>2 = ^(X fe - 7 X*-i-p) 2 . (u 



k=l k=l 

In particular, those estimators of (b, a) are given by 

% n = _ log ELi ELi Xk-n ELi Xk-iX* (19) 

and 



1 _ e -fen; 



&n- (1.10) 



Following Wei and Winnicki (1990), we also consider the WCLSEs of (7, p) and {b, a) by 
minimizing the weighted sum 

4 _^[X*-7(**-i + 1)-(p-7)] s qir 



n 2 « 

f-f X fc _! + 1 f-f 



+ 1 



The reason of considering the above quantity is it does not fluctuate too much even when 
the values of the samples X k , k — 0, 1, • • • , n are large. The resulting WCLSEs of {b, a) 
are given by 

v 2^ k =i^k 2^ k =i x k _ 1+ i ,l ^=ix H +i , v 

6„ = - log , rTTv^ 1 2~~ L12 

E fc =i(^fc-i + 1) Efc=i ^rr+T - 72 
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and 



E n Y — P~ h n Y 

a n — - u n . 

n(l — e ° n ) 

The main purpose of this paper is to study the asymptotic properties of the CLSEs 
and the WCLSEs given above. We show that the estimators are consistent and obey some 
central limit theorems. In particular, for 1 < a < 2 we prove that n^-^^ibn-b, a n — a) 
converges to an a-stable random vector as n — > oo. For 1 < a < (1 + v5)/2, we show 
that n {a - l)/a \b n - b, a n — a) converges to a nontrivial limit as n — > oo. A combination 
of this with the result of Overbeck and Ryden (1997) for a = 2 only gives a partial 
characterization of the asymptotic behavior of the CLSEs. The characterization of the 
exact asymptotics of the CLSEs for (1 + \/5)/2 < a < 2 is left as an open problem. 

The proofs of our limit theorems are very different from and much harder than the 
Gaussian case. The key of the approach is to establish the convergence of some point 
processes and partial sums associated with a stationary realization of the SCIR-model. 
The techniques in this subject have been developed extensively by Basrak and Segers 
(2009), Davis and Hsing (1995), Davis and Mikosch (1998) among others. We also make 
use of the results of Hult and Lindskog (2007) on the extremal behavior of Levy stochastic 
integrals. The proofs depend heavily on the construction and characterization of CBI- 
processes in terms of stochastic equations of given in Dawson and Li (2006, 2012), Fu and 
Li (2010) and Li and Ma (2008). 

We finally propose an estimator of the volatility coefficient a based on high frequency 
observations at times {0, 1/n, - ■ • ,(n — l)/n, 1}. Suppose that the parameter a is known. 
Given constants p G (0,a) and 5 G (0,min{l — 1/a, 1/a 2 }), let 



1 / I Xk/n — X(k-l)/ n 



fc=l A (fc-l)/n 

We prove that a n is a weakly consistent estimator for a. 

The paper is organized as follows. In Section 2, we prove the exponential ergodicity 
of some subcritical CBI-processes, which implies the strong mixing property of the SCIR- 
model. Section 3 is devoted to the regular variation properties of some random sequences 
defined from the model. The limit theorems of random point processes and partial sums 
are established in Section 4. Based on those theorems, the asymptotic properties of the 
estimators are proved in Section 5. 

Notation. Let N = {0, 1, 2, • • • } and Z = {0, ±1, ±2, • • • }. Let K = (-00,00), 
R = [-00, 00] and Kg = R d \ {0}, where = (0, 0, ■ ■ • , 0). Let C + (M$) be the collection 
of positive continuous functions on Rjj with compact support. Let M(Rq) be the class of 
Radon point measures on Rq furnished with the topology of vague convergence. We use C 
with or without subscripts to denote positive constants whose values are not important. 
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2 CBI-processes and ergodicity 



In this section, we prove some simple properties of CBI-processes. In particular, we prove 
a subcritical CBI-process is exponentially ergodic and strongly mixing. The results are 
useful in the study of the asymptotics of the estimators. We start with an important 
special case of those processes. Let a > and b be constants and (u A u 2 )m(du) a finite 
measure on (0, oo). For z > set 



oo 



(z) = bz + -a 2 z 2 + / (e~ zu - 



o 



A Markov process with state space M + := [0, oo) is called a continuous-state branching 
process (CB-process) with branching mechanism <p if it has transition semigroup (Qt)t>o 
given by 



oo 

e~ > *Q t {x,dy)=e- m *W, (2.1) 



o 

where t \-> v t (X) is the unique positive solution of 

JU t (A) = -</>(y t (\)), v (X) = X. (2.2) 

The CB-process is called critical, subcritical or supercritical as b = 0, b > or b < 0, 
respectively. From (12. 2p we obtain the following semigroup property: 

v r+t (X) = v r (v t (X)), r,t,X> 0. 

Taking the derivatives of both sides of ( \2.2h one can see u t := (d/d\)v t (0) solves the 
equation (d/dt)u t = —bu t , and so u t = e~ bt for t > 0. Then differentiating both sides of 
flZZE]) gives 

^"OO 

yQ t {x,dy) = xe- b \ t,x>0. (2.3) 



n 

By Jensen's inequality, we have v t (X) < Xe~ bt for t, X > 0. 

It is easy to see that (Qt)t>o is a Feller transition semigroup, so it has a Hunt realiza- 
tion. Let X = (f2, £f , £ff, X t , Q x ) be a Hunt realization of the CB-process. The hitting time 
To = inf{t > : X t = 0} is called the extinction time of X. It follows from Theorem 3.5 
of Li (2011) that for t > the limit v t =f liniA->oo v t(X) exists in (0, oo], and 

Qx(r <t) = Q x (X t = 0) = exp{-x^}. (2.4) 

By Theorem 3.8 of Li (2011), we have v t < oo for all t > if and only if the following 
condition holds: 

Condition 2.1 There is some constant 6 > such that <p(z) > for z > 9 and 

(biz^dz < oo. 
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Proposition 2.2 (Li, 2011, p. 61) If Condition \2. 1\ holds, then v t = l t (0,oo) is the mini- 
mal solution of 

^-v t = -<j>(v t ), v = oo. (2.5) 

Let t i — y v t (X) be defined by (12. 2p . A Markov process with state space M + is called a 
CBI-process with branching mechanism <fi and immigration rate a > if it has transition 
semigroup (P t )t>o given by 

-xv t (\)-aj v s (\)dsy (2.6) 
By differentiating both sides of (12. 6p we obtain 

rOO rt 

/ yP t (x, dy) = xe~ ht + a e~ hs ds = xe~ ht + ab~ l (l - e~ bt ), (2.7) 
Jo Jo 

where b~ l (l — e~ bt ) = t when b = by convention. 

A realization of the CBI-process can be constructed as the strong solution to a stochas- 
tic integral equation. Let W(ds, du) be a time-space Gaussian white noise on (0, oo) 2 with 
intensity dsdu and N\{ds,dz,du) a Poisson random measure on (0, oo) 3 with intensity 
dsm(dz)du. Let Ni(ds,dz,du) = N^ds^dz^du) — dsm(dz)du denote the compensated 
measure. Then for each x > there is a pathwise unique positive strong solution to the 
following stochastic equation: 



rt rt rY s -{,x) 

Y t (x) = x + / (a - bY s (x))ds + a / W(ds, du) 

Jo Jo Jo 

_(*) 

zN^ds^dz^du). (2.8) 



JO JO 



The solution {Y t (x),t > 0} is a CBI-process with branching mechanism and immigration 
rate a. See Theorem 3.1 of Dawson and Li (2012) or Theorem 2.1 of Li and Ma (2008). 
A slightly different formulation of the process was given in Dawson and Li (2006). 

Proposition 2.3 Suppose that Condition \2.1\ holds. For x,y > let T x ^ y := inf{t > : 
Y t (x) — Y t (y) = 0}. Then we have P{T x y < oo} = 1 and 

P{T x , y <t} = exp{-|x - y\v t }, t > 0. (2.9) 

Moreover, we have Y t (x) = Y t (y) for all t > T x y . 

Proof. It suffices to consider the case of y > x > 0. By Theorem 3.2 of Dawson and Li 
(2012), we have P{Y t (x) > Y t (y) > for all t > 0} = 1 and {Y t (x) - Y t (y) : t > 0} is a 
CB-process with branching mechanism 0; see also Remark 2.1 (iv) of Li and Ma (2008). 
Then (12.91) follows from (12.41) . The pathwise uniqueness of (12.81) implies that Y t (x) = Y t (y) 
for all t > T XiV . By Corollary 3.9 of Li (2011) we have P{T XtV < oo} = 1. □ 
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The above proposition provides a successful coupling of the CBI-processes. This has 
many important implications. We refer the reader to Chen (2004) for systematical study of 
coupling methods and their applications in the theory of Markov processes. In particular, 
we shall use the above coupling to prove the strong Feller property and exponential 
ergodicity of the CBI-process following Chen (2004, p. 37). Write / G b&(R + ) if / is a 
bounded measurable function on R + . 

Theorem 2.4 Under Condition \2.1\ the transition semigroup (Pt)t>o given by $2. 6}) has 
the strong Feller property. Moreover, for any t > 0, x, y > and f G b=^(R + ) we have 

\PJ(x) - PJ(y)\ < 2U/IU1 - e-*l*-»l), (2.10) 

where \\f\\oo — su Pa; 1/0*0 1 denotes the supremum norm. 

Proof. It suffices to prove (12.101) . By Proposition 12.31 we have 

\PJ(x) - PJ(y)\ = \E[f(Y t (x)) - f(Y t (y))]\ < 2||/|| 00 P(T., 8 > t), 
which together with (12.91) implies (12.10p . □ 



Theorem 2.5 Suppose that b > 0. Then the transition semigroup (Pt)t>o has a unique 
stationary distribution \i, which is given by 

Lp(\) = J e~ Xx ^{dx) =exp{ -ajf z^z)' 1 dz} , A > 0. (2.11) 



Moreover, we have 



d 



d\ 



x^dx) = —L,(X) _ = -. (2.12) 



a=o+ b 



Proof. By Theorem 3.20 of Li (2011) and its proof given there, for any x > we have 
/i = limt-^oo Pt(x, ■) by the weak convergence. Then (Pt)t>o h as the unique stationary 
distribution \i. By differentiating (12. lip we obtain ( 12. 12ft . □ 

Theorem 2.6 Suppose that b > and Condition \2.1\ is satisfied. Then the transition 
semigroup (Pt)t>o is exponentially ergodic. More precisely, for any x > and t > 1 we 
have 

\\P t (x, - MOIU < 2(1 - expj-^xe- 6 ^- 1 )}) + 2v ia b~ l e"^, (2.13) 
where \i is given by 112.11}) and \\ ■ \\ V3X denotes the total variation norm. 

Proof. We only need to prove (I2.13p . In view of (12.111) . we have 

< 1 - L M (A) < a [ z(j){zY l dz < a6 _1 A. (2.14) 
Jo 
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Then, for any / G hi 



with 11/11 < 1, we can use Theorem 12.41 and (12.141) to see 



\Ptf(x) - Kf)\ < I \Ptf(x) - PJ(y)\n(dy) < 2 

'o Jo 



-v t \x-y\ 



)fi(dy) 



<2 / (1 
Jo 

< 2(1 -e 



~'Mdy) + 2 



-vty 



)fj.{dy) 



-**) + 2 [1 - L,{v t )} < 2(1 - e"**) + 2ar 1 ^. 
By Proposition I2.2[ for t > 1 we have w 4 = Ut_i(iJi) and so v t < e~ b ^^Vi. Then we obtain 



Under the conditions of Theorem I2.6[ for any finite set {t\ < t 2 < ■ ■ ■ < t n } C Kwe 
can define the probability measure ^ t 2 ... t n on ^+ by 



, dXnj 



/J't 1 ,t 2 ,-,t n (dx 1 ,dx 2 , 

= fi(dx 1 )P t2 - tl (x 1 ,dx 2 ) ■ ■ ■P tn -t n _ 1 (x n -i,dx n ). 



(2.15) 



It is easy to see that {^t lt t 2 ,- ,**» ■ *i < *2 < ■ • • < t n G M} is a consistent family. By 
Kolmogorov's theorem, there is a stochastic process {Y t : t G M.} with finite-dimensional 
distributions given by (I2.15p . This process is a (strictly) stationary Markov process with 
one-dimensional marginal distribution /i and transition semigroup (Pt)t>o- Since (Pt)t>o 
is a Feller semigroup, the process {Y t : £ G M} has a cadlag modification. 

Theorem 2.7 Let {Yf : t G M} fre a Markov process with finite- dimensional distributions 
given by A2.15\) . Then it is strongly mixing with geometric rate, that is, as t — >• oo, 

7r t := sup sup \P{AnB)-P(A)P(B)\ 

Aea{Y s ,s<0} Bea{Y 3 ,s>t} 

decays to zero exponentially. 

Proof. It follows from (l2~IHj) and (12141 that, for t > 1, 

/■oo 

||P t (z, •) - < 2 / (1 - e-^-^Mda;) 

+ 2u 1 afe- 1 e- 6( *- 1) 

< 2[1 - ^(vie-^*-^)] + 2v 1 ab- 1 e~ b{t - 1) 

< 4D 1 ar 1 e- f ' (t - 1) . 

Then {Y t : t G R} is strongly mixing with geometric rate; see, e.g., Mikosch and Strau- 
mann (2006, p.516) or Bradley (2005, p. 112). □ 



Lemma 2.8 Suppose that {Z t } is an a-stable Levy process with 1 < a < 2 and {y(t)} is 
a predictable process satisfying, a.s., 

-T 

\y(t)\ a dt < oo, T>0. 



Then for any < r < a, there exists a constant C = C(r, a) > such that 



E 



sup 

L t<T 



y{s)dZ< 



< CE 



T \ r/a 

\y(t)\ a dt 
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Proof. When {y(t)} is a positive process, the result follows from a result proved in Long 
and Qian (2011). Then, in the general case, we have 



E 



sup 

t<T 



y(s)dZ s 



r~ 








r 




r~ 




< CiE 


sup 


/ y+(s)dz s 


+ sup 


/ y~(s)dZ s 








L t<T 


Jo 


t<T 


Jo 





< C 2 E 

< 2C 2 E 



T 



y+{t) a dt 

T \y(t)\ a dt 



r/a 



+ 



y-{t) a dt 



r/a 



r/a 



where y + and y_ denote respectively the positive and negative parts of y. 



□ 



Now let us consider a filtered probability space (Q, J&" t , P) satisfying the usual 
hypotheses. Let {Z t : t > 0} be a spectrally positive a-stable Levy process. For a = 2 we 
understand the process as a standard Brownian motion; and for 1 < a < 2 we assume it 
is a stable process with Levy measure v a {dz) given by ( 11 .3(1 . By Theorem 6.2 of Fu and Li 
(2010), for any initial value Xq, which is a positive J^cr m easurable random variable, there 
is a unique positive strong solution {X t : t > 0} to (1 1 . 2 (1 . The existence and uniqueness 
of this solution also follows from Corollary 6.3 of Fu and Li (2010) by a time change. Let 
/ be a bounded continuous function on K with bounded continuous derivatives up to the 
second order. For a = 2, we can use Ito's formula to see that 

f(X t ) = f(X r ) + y* Lf(X s )ds + M t (f), t > r, (2.16) 

where {M t (f) : t > r} is a martingale with respect to the filtration {,^t)t>r and 

Lf(x) = (a — bx)f'(x) + —xf"(x), x > 0. 

When 1 < a < 2, by the Levy-Ito representation of we can rewrite (11.2)1 into the 

integral form: 



X t = X r + / (a- bX s )ds + tr 



t poo 



r JO 



XlL a zN{ds,dz), t>r, 



(2.17) 



where N(ds, dz) is a compensated Poisson random measure on (0, oo) 2 with intensity 
dsu a (dz). By Ito's formula one can see that (12. 16(1 still holds for 1 < a < 2 with the 
operator L defined by 



Lf(x) = (a- bx)f'(x) + 



a 



aT(— a) 



\f{x + y)-f{x)-yf{x)\ 



y 



dy 

a+l ' 



By Theorem 9.30 of Li (2011), for any 1 < a < 2 we can identify the SCIR-model as a 
subcritical CBI-process with immigration rate a and branching mechanism 



a 



0(A) = b\ + — \ c 
a 



A > 0. 
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It follows from Theorem 12.51 that the SCIR-model has the unique stationary distribution 
\x with Laplace transform given by 



L »( X ) = J q e- Xx fx(dx) = exp { - J 



aadz 



ab + a a z a ~ 1 



A > 0. 



(2.18) 



Let P x denote the law of the SCIR-model {X t : t > 0} defined by flT2) with X = x > 
and let denote the corresponding expectation. 

Proposition 2.9 Suppose that < a < 2. Then for any < (3 < a, there is a constant 
C > so that, for t, T > 0, 



E, 



(\f e ~ Kt ~ S)XX s- dZ ^\) 



x /3/a e -0bt/a 



) 



and 



E x ( sup 

\ 0<t<T 



e -b(t-s) x l/a dz 



Proof. Using Lemma [2 .8} we have 



E, 



* \P/a 

e abs XJs 



o 



and 



E x ( sup 

\ 0<t<T 



-Kts) x l/a dz 



< CxE, 



T , \ (8/a 

e abs XJs 



o 



By Holder's inequality and (12.71) it is easy to see 



E, 



o 



e abs X s ds) 



< 



< 



Jo 



1 u -\P/ a 
E x (e abs X s )ds 



e abs (xe~ bs + a&" 1 )ds 



/3/a 



Then we have the desired inequalities. 



□ 



Proposition 2.10 Suppose that < a < 2. T/ien /or any < f3 < a, there is a constant 
C > and a locally bounded function T i— > C(T) > so £/ia£, /or £, T > 0, 



and 



E x ( sup Xf) < C(T)(l + x* 
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Proof. Using (ll.5p with r = and an elementary inequality, we have 



and 



-6(t-s) T^ 1 /" 



E,( sup Xf ) < d 

■ 0<t<T 



+ a^b~P + er^E x . ( sup 

0<t<T 



Then the results follow by Proposition 12.91 



' e- h ^X]L a dZ, * 



□ 



3 Regular variations 

In this section, we study the regular variation property of some random sequences asso- 
ciated with the SCIR-model. We shall make use of the stochastic equations ( 11.21) and 
(11.51) as well as the results of Hult and Lindskog (2007) on the extremal behavior of Levy 
stochastic integrals. The reader may also refer to Samorodnitsky and Grigoriu (2003) for 
some results on the tail behavior of solutions to certain stochastic differential equations 
driven by Levy processes. For the convenience of the reader, we first recall some concepts 
and properties of regularly varying sequences. Let "| • |" be any norm on M. d . 

Definition 3.1 A d-dimensional random vector X is said to be regularly varying if there 
exists a Radon measure t] on M d , finite on sets of the form {x G M d : \x\ > r}, and a 
sequence {a n } satisfying a n —> oo such that, as n — > oo, 



nP(a" x X G 



(3.1) 



The above sequential form of the condition is the same as saying there exists a (necessarily 
regularly varying) function t \-> g(t) such that, as n — > oo, 

It is known that the condition implies the existence of a constant a > such that t] (r A) = 
r~ a r](A) for all r > and all A G 3S(M. d ) bounded away from 0, where rA = {rx : x G A}. 
In this case, we say X is regularly varying with index a > 0. 

Proposition 3.2 Let £ be a positive regularly varying random variable with index a > 0. 
Then we have: 



[i) If a > 1, then as x — oo, 



E(£l{c>x}) 



a 



a — 1 



xP(£ > x); 



(ii) IfO < a < 1, then as x — >■ oo, 

E(ei {? <x } )~ r ^xP(e>x). 



ii 



Proof. The results are immediate consequences of Karamata's theorem; see Resnick (2007, 
p. 25). Let G(x) — P(£ > x). By integration by parts, 

/*oo /*oo 

E(ei {€ >, } ) = - / ydG(y) = xG(x) + / G(y)dy. 

J X J X 

The second term on the right hand side is equivalent to (a — l)~ 1 xG(x). Then (i) follows. 
The proof of (ii) is similar; see also Resnick (2007, p. 36). □ 

There are several equivalent forms of the regular variation property; see, e.g., Resnick 
(1986, p. 69). One of them is given in the next theorem. Let S d_1 = {x G M d : \x\ = 1} be 
the unit sphere. 

Theorem 3.3 A d- dimensional random vector X is regularly varying with index a > 
if and only if there exists a probability measure v on S d_1 such that for every r > 0, as 
t — > oo, 



P(|X| > rt, \x\- l x e 
P(|X|>t) 



r~ a u(-) 



Definition 3.4 A sequence of random variables {X^ : k G Z} in R d is called jointly 
regularly varying if all the vectors of the form (Xi, ■ • • , X;) are regularly varying. 

In the sequel, we shall also use the norm ||x|| := maxj|xj| for x = (x±, ■ ■ ■ ,Xd) G 
~R d . By Corollary 3.2 in Basrak and Segers (2009) we have the following equivalent 
characterization of the jointly regular variation property: 

Theorem 3.5 Let {X& : k G Z} be a stationary sequence in M. d and assume x i— > 
P(||Xo|| > x) is regularly varying with index —a for some a > 0. Then {X^ : k G Z} is 
jointly regularly varying with index a if and only if there is a sequence {@k '■ k G N} such 
that for every k G N, as x — > oo, 

PfllXoir^Xo,--- ,X fc ) G -HlXoll >z) -^»P ((©(,,••• ,0 fc ) G •). 

The concept of regular variations can also be defined for continuous time stochastic 
processes. Let T > and let D d [0, T] := ©([0, T], M. d ) be the space of all Revalued cadlag 
functions on [0, T] equipped with Skorokhod topology. Let 



S d [0,T] = \ye3 d [0,T]: sup ||y t || = l). 

0<t<T 



Definition 3.6 A stochastic process Y = {Y t : < t < T} with sample path in D d [0, T] 
is said to be regularly varying if there exist a measure Q on D d [0,T], finite on sets 
bounded away from 0, and a sequence {a n } satisfying a n — > oo such that for any set 
B G ^(D d [0, T]) bounded away from with Q{dB) = 0, as n -> oo, 

nP(a^ 1 Y eB)^ Q(B), 

The above property implies there is a constant a > such that Q(uB) = u~ a Q(B) for 
all u > and all B G £$(B d [0, T]) bounded away from 0. In this situation, we say Y is 
regularly varying with index a > 
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The convergence in the above definition can be formulated for general boundedly finite 

measures on 6g[0,T] = (0, oo] x § d [0,T]. We shall denote the convergence by "— — >■". The 
reader may refer to Hult and Lindskog (2005) for more details. 

Remark 3.7 Let < a < 2 and let {Z t : t > 0} be a one-dimensional a-stable process 
with Levy measure v{dz). It follows from Lemma 2.1 of Hult and Lindskog (2007) that, 
as n ^ oo, 

nP{n^ a Z t G ■) —> tu(-). 

Remark 3.8 Let < a < 2 and T > 1. Suppose that {Z t : < t < T} is a one- 
dimensional Levy process such that X — Z\ satisfies 113.1}) with 77(21, 00) = cz~ a for some 
c > 0. Let {Y t : < t < T} be a positive predictable caglad process satisfying sup 0<t<T Y t > 

a.s. and E[sup 0<t<T l^ a+<5 ] < 00 for some 5 > 0. By Theorem 3.4 and Example 3.1 in 
Hult and Lindskog (2007), for any z > and < t < T we have, as n — > oo ; 

nP[a7 n l J Y s _dZ s > z^j -> V (z, 00) J E{Y s a )ds = cz~ a J E{Y s a )ds. 

Remark 3.9 Suppose that {Y^} is a stationary sequence of regularly varying random 
vectors. Let {a n } be taken such that nP(|Yo| > a n ) — > 1 as n — >■ 00. By Lemma 2.3.9 
of Basrak (2000), the strong mixing condition implies the mixing condition £/(a n ), i.e., 
there exists a sequence of positive integers r n such that r n — > oo ; l n = [n/r n ] — > 00 as 
n — > 00 and 

n r„ , 

Eexp { - /(Yfc/fln)} - (Eexp { - ^ /(Y fe /a„)}) " 
fc=i fc=i 

/or every / G Cq"(Ro)- / ac ^ «/ {Y&} strongly mixing with geometric rate, we can 
choose r n = [n s ] for any < 5 < 1/ see Remark 2.3.10 of Basrak (2000). 

Some regularly varying distributions can be found in the SCIR-model. Recall that 
P x denotes the law of the SCIR-model {X t : t > 0} defined by (JO} with X = x > 
and E a denotes the corresponding expectation. In the sequel of this section, we assume 

1 < a < 2. 

Proposition 3.10 For any x > we have, as u — > 00, 

P x (X t > u) ~ ° - [g a (t)+p Q (t)x]n" a , 
al (— a) 



where 



' -bt „-abP a 



1 



P«(*) = I7-^T[e" M - e' ab % q a (t) — y —7(1 — e^) -p a (t) . (3.2 
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Proof. In view of ( II. 5p . the extremal behavior of X t is determined by a stochastic integral. 
Then, using Remark 13. 8| we have, as u — > oo, 

P x {X t >u)~ PxfaJ e"' ( '" s) 0^dZ s > u) 

~ o^P^ > u) [ e- abit - s) E x (X s )ds. 
Jo 

Based on (12 .7p . it is easy to compute 

e -<*Kt-s) Xsds \ =q a (t)+p a (t)x. (3.3) 



-t 

E, ' 



By Remark 13.71 we have ~P x {Z t > u) ~ t/acT(— a)u a . Then the desired result follows. □ 
Proposition 3.11 For any K > 0, we have 



lim sup 



u a P s (Xi > u) — -(q a +p a x) 

ai (—a) 



and 

0. 



lim sup 

u ->°° xe[o,K] 



u a ^(X^x^y) - — — Aq a +p a x) 

[a — 1)1 —a 



where p a = p a (l) and q a = q a (l) are defined by 113. ty) . 

Proof. Let {X t (x)} be the SCIR-model defined by ( 11.21) with initial value Xq = x. 
By Theorem 5.5 of Fu and Li (2010), the random function x i— > X t (x) is increasing, 
so x H- P x (X t > u) is increasing for any t, u > 0. Then the first convergence holds 
by Proposition 13.101 and Dini's theorem. The second convergence follows similarly by 
Proposition 13.21 □ 

Let us consider a stationary cadlag realization {X t : t G M} of the SCIR-model with 
one- dimensional marginal distribution \i given by (12.181) . By a modification of the argu- 
ments in the proofs of Theorems 9.31 and 9.32 Li (2011) one can see that, on an extension 
of the probability space, there is a compensated Poisson random measure N(ds, dz) on 
R x (0, oo) with intensity dsu a (dz) so that (I2.17P is satisfied for all t > r G R. For any 
integer k G Z let 

I k = s k (l, (1 + Xfc.!)" 1 ), H fc = X^Xl^Bk), (3.4) 

where 

pk poo 

e k = a / e- b{k - s) X l J_ a zN{ds,dz). (3.5) 

Jk-l Jo 

It is easy to see that the above sequences are stationary. We are going to prove that the 
sequences {X k }, {l k } and {H fc } are jointly regularly varying. 
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Proposition 3.12 Let \i be the stationary distribution of the SCIR-model given by A2.18\) . 
For any t > we have, as x — > oo, 



fi(x, oo) 



aa 



a 2 b 2 



r(i-a) _ V 



a 3 b 2 



T(-a)- 1 x- a . 



Consequently, for any < r < a we have 

/"OO /"OO 

/ x r [i(dx) = / ^{y 1 ^ ', oo)dy < oo. 
Jo Jo 

Proof. The tail behavior of X 4 is closed related with the asymptotics of its Laplace 
transform. By (I2.18p . as A -»■ 0, 



r ... , a /" A a fate , 
L M A=1--/ — - + 0(A 
b J ab + a a z a 1 



a a 
1-tA + t 



— + ^ 



1 - 



a, a f x ba a z a ~ l dz ^ /l2 , 
6 A + FX at+ff ^- +0 < A ' 



A 



b ab 2 J ab 2 J ab + a a 2; c 

1 - A / + ^A Q - 0(A 2a_1 ) + 0(A 2 ) 



a 



where we have used (12.121) for the last equality. Then the result follows by Theorem 8.1.6 
of Bingham et al. (1987). □ 



Lemma 3.13 Let N(ds,dz) be the compensated Poisson random measure in ( 2.11 ) and 
let 



-t pi 
z(t) = al I e 



JO 



b ^xlL a zN(ds,dz). 



(3.6) 



Then for any 1 < r < a 2 and any T > 0, we have 



E 



sup \z(t)\ 7 

0<t<T 



< oo. 



Proof. We follow an idea in the proof of Lemma 5.5 of Hult and Lindskog (2007). Let 
q = 1/(1 — p~ x ) for any p G (l,a 2 /r). By the Burkholder-Davis-Gundy inequality and 
Holder's inequality, we have 



e bs xlL a zN(ds, dz) 
e 2hs X 2 J_ a z 2 N(ds ) dz) V ' 



E 


sup \z(t) \ r 


< a r E 


sup 


fj 




-0<t<T 




Lo<t<T 


Jo Jo 



< CiE 

< C(T)E 



7 

° T 1 

sup X r J a ( [ [ z 2 N(ds,dz) 



r/2 
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< C(T)E 1/p sup X r / /a E l/q 

-0<s<T 



T pi 



JO 



z 2 N(ds, dz] 



rg/2-. 



The first expectation on the right-hand side is finite by Proposition 12. 101 By Theorem 34 
of Protter (2005, p. 25), the second expectation is also finite. □ 

Lemma 3.14 Suppose that {A n } C J^o is a sequence of events so that P(A n ) — > as 
n — > oo. Then for any x > and T > we have 



> a n x ) = 0. 



limnP(A n , sup / e~ b{t - s) X l J_ a dZ s 

n ^°° ^ 0<t<T Jo 

Proof. The right-hand side of (13. 7p is bound above by J\ + J<i + J3, where 



(3.7) 



Ji = nP ( sup 

■ 0<t<T 



Jo 

-t 



e- b ^XlL a zN(ds, dz) > a n x/3 



and 



J 2 = nP( sup / e- b{t - s) X] /a ds ! zu a (dz) >a n x/'i\ 
J 3 := nP(A n , sup / f e" 6( *- s) 'X^ Q zN{ds,dz) >a n x/3) 

\ 0<t<TJo Jl ' 



Let z(t) be defined by (13.61) . Then for any 1 < r < a 2 we have 

sup |z(t)| r 



4<-S=-E 



0<<<T 



X 1/Q ds 



where the two expectations are finite by Proposition 12.101 and Lemma 13.131 Then J\ — > 
and J 2 — > as n — > 00. By introducing the Levy process 



t /"OO 



f(t) : = y y zN(ds,dz), t > 0, 

for any > 1 we have 

J 3 < nP(A n ,i^/ a £(T) > a n x/6) + nP (J o Xl-l{x.->K}d£(s) > a n x/Q). 

By Remark 13.71 and the property of independent increments of {£(£)} it follows that 
lim nP(A ni K l/a i{T) > a n x/6) = lim nP(A n )P(K 1/a £(T) > a n x/Q) = 0. 

By Remarks 13.71 and 13.81 

UranP^y Xl[ a l {Xs _ >K} d£(s) > a n x/6^) = C 3 x~ a E[X s l {Xs>K} ds, 
which tends to zero as A' -)• 00. Then we have the desired result. 



□ 
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Theorem 3.15 The sequence {X^} is jointly regular varying with index a. More pre- 
cisely, as x y oo we have 



P(X > x 



a(T TV \-l 

^ — r r(— a) x 



a 3 b 2 ' 



and, for any integer k > 1, 

p(x - 1 (x ,--- ,x k )e-\x >x) 



. t e-bk)\ 



(3.8) 



(3.9) 



Proof. By Proposition 13.121 we have the asymptotics (13. 8p . It suffices to show (13. 9p 
holds when x —t oo along the sequence a n := n 1//a . Let X = (X ,Xl,--- ,Xk) and 
X = (Xo,X e~ b , ■ ■ ■ ,X e~ bk ). Let z{t) be defined by (13. 6p . For any 5 > 0, we can use 
ffTBD to see 



P(||X-X|| > a n 5|X > a n ) 



< P ( ak + max 

i<j<fc 



e -bV-') x V a dZ, 
'■.I 



> a n b~ 



< CnP\Xr, > a n , max 



o 



X > a n 
> a n 5 — ak 

> (a n 5 — ak)/k) . 



By Lemma T3. 141 it is easy to see the right-hand tends to zero as n — > oo. By Theorem 13.31 
we conclude that, as n — > oo, 



P(X/a„ G -\X > a n ) ~ P(X/a n G -|X > a T 



a 



dz 



-{z(l,e-»,- ,e-»)e-} 7^1- 



Then (13.91) follows by the continuous mapping theorem. By Theorem 13.51 the sequence 
{Xfc} is jointly regular varying with index a. □ 



Theorem 3.16 The sequence {I&} defined by IjS.fy is jointly regular varying with index 
a. More precisely, as x — >■ oo we have 

and, /or any integer k > 1, 

Pfllliir^Ii,--. ,h) e -llllill >z) 

^> E(G)- 1 E [G; ((1, (1 + Xo)" 1 ), 0, • ■ ■ , 0) G •] , (3.11) 

G = a a [ e~ ab ^X t dt. (3.12) 



17 



Proof. Clearly, it suffices to show that ( 13. lUj) and ( 13. lip hold when x — > oo along the 
sequence a n := n l / a . For < s < k, let Z s = (Z S ,Z S ) and <fr s = (0i(s), (f> 2 (s)), where 

Ms) = J2 1 u-ids)e- b{j - s) °x 1 s- a i Ms) = Eiy-mOO— 

j=i i=i 



We consider the process 



i 2 (s)dZ s ). 



Choose some 6 G (a, a 2 ). By Proposition 13.121 we have E[Xq^*] < oo. Then Proposi- 
tion 12.101 implies that 



E 



sup ||$« 



0<s<k 



< E 



sup a 

0<s<k 



8 y-8/a 



< a d C{k)[l + E{X S /a )} < oo. 



By Theorem 3.4 in Hult and Lindskog (2007) and Remark 13.71 as n — > oo, 

nP(a- l (® ■ Z) G •) A Q(-) := kE[u a {x G M+ : x$ T l[ T ,k] G ■}], (3.13) 

where ^ Q is defined by (ll.3p and r is uniformly distributed on [0, k] and independent of 
In view of (13.131) . we have 

nP(a-\$ ■ Z)i G ■) Q({y G O 2 [0, fc] : yi G ■}). 

By Definition 13.11 it follows that 

Q({y G D 2 [0, fc] : ||yi|| > r}) = r~ Q Q({y G D 2 [0, fc] : ||yi|| > 1}). 

Now define the functions h ,hi,h 2 : D 2 [0, fc] — > M 2fc by 

My) = (yi, y2 - yi, • • • , y& - y^-i), 
My) = ^iiyii^i}' ^(y) = h o(y)h(y)- 

Let Disc(/ij) be the set of discontinuities of hi (i = 0,1,2). By (13.131) it is easy to see 
that Q(Disc(feo)) = Q(Disc(/ii)) = 0, so Q(Disc(/i 2 )) = 0. Moreover, for any B G M(R 2k ) 
bounded away from the set h 2 1 (B) G £g(B 2 [0,k}) is bounded away from 0. Applying 
the continuous mapping theorem, we obtain as n — > oo, 

nPdllxll > a n , a-^hJz, ■ ■ ■ , I k ) G •) Q o h^(-) 

on M 2fc \ {0}, where 

Qoh 2 \-) = m[v a {x G M+ : ||x* T l [Tifc] (l)|| > 1, 

s(* T l [T|fc] (l),0,-- - ,0) G ■} 



kE 



v a \x G 



: I la;*- 1 1 > 1, r < 1, 
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Let E = {(y 1; 



x(* T l [Tjfc] (l),0,-.. ,0) G •} 
= kE u a {x G R+ : x\\$ T \\ > l,x($ r ,0, • • • , 0) G ■}l [0i i](r) 

, yjt) G R 2k : ||yi|| > 0}. Define the injection / : E -+ (0, oo) X E by 

/(yi,--- ,yfc) = (||yi||,yi/||yi||,-- - ,yfc/||yi||)- 



Then we have as n — > oo, 

nPdlixll > a n , Hiiir^li, -,i*)e-)^Qo hr 2 x o r l (-) 

on where 

v a {xeR + :x\\& T \\ > l,(* T /||* r ||,0,--- ,0) G ■}l [0 ,i](r) 
v a {xeR + ■■x\\®t\\ > l}l [0 ,i](r),(* T /||* T ||,0,.-- ,0) G ■ 



kE 
kE 



k r roc 
E 



aT(—a 
k 



dz 



a+1 



l [0 ,i](r),((l,(l+X )- 1 ),0,-- - ,0) G 



-EOl^flp^T), ((1, (1 + Xo)" 1 ), 0, • ■ ■ , 0) G ■] 

^— / E[||$ s |r,((l,(l + X )^),0,---,0)G-]rf S 
I a ) Jo 



a 2 T( 



a 2 T(—a) 
In particular, as n — > oo we have 



E[G, ((1,(1 + Xq)- 1 )^,--- ,0) G •]. 



nPflllJ > a r 



E(G) 

a 2 r(— a) 



By (gUD , we have E(X t ) = a/b. It follows that 
E(G)= le~ 



acr" 



(1 -e" a6 ). 



Then we have f 1 3 . X U j) and (13.111) . By Theorem 13.51 the sequence {1^} is jointly regular 
varying with index a. □ 



Remark 3.17 For k = 1, 2, • • • define 



V k = a e- b ^e~ b ^~ k+1 ^ a dZ, 



(3.14) 



fc-i 

Then the sequence {14} is i.i.d. with the same distribution as 

'e~ b -e~ 



CT 



-b _ „-ab^ i/ a 



(a - 1)6 

which is regularly varying with index a. 



Zi, 
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Lemma 3.18 Let V k be defined by \3.1J$ and let H k = ^"^""(l, V&). Then for any 
< r < a 3 / (a 2 + 1), we have 

E[||H fc -H fc || r ] <oo. 
Proof. Since < r < a 3 (a 2 + < a, by Lemma [2.81 and Holder's inequality, 

E[||H fc -H fe |r] =E 



k 
fc-1 



< CE 



ae- b{k - s) X k ^xlL a - X 1 k L a 1 e- b(s - k+1)/a )dZ ; 

-ab(k-s) X *\ Xa _ X fc „ ie " 6(s " fc+1) |rfs V 



< 
< 



CE{xi_ 1 [E Xk _^^ \X S - X e- bs \ds)] T/a } 
CE{x r k _ 1 [j^ E Xk ^ (\X S - X e- bs \)ds] r/a } 



< CE{Xl_, 

< ceIxu 



E 



e -b(s-u) x l/* dZu 
e -Hs-u) X l/° dZu 



ds 
ds 



r/ot 
1/a 



< ce[xu{i + xHX 

which is finite by Proposition 12.101 



□ 



Theorem 3.19 The sequence {H^} defined by jjj3.4\) is jointly regular varying with index 



a 



''/(a + 1). Let V\ be defined by ^3.14 )- Then we have, asx-> oo, 



PdlHj >x) -EflViyxl 02 /^] 



^_ x -a 2 /(a+l) 



a 3 b 2 T(— a) 



and, for any integer k > 1, 



where 



P(||H 1 ||- 1 (H 1 ,...,H fc )G.|||H 1 ||>a:) 

w E[iv|v 1 | aa /( a + 1 ) ; (e 1> ---,e fc )€-] 

~* E[lV|^i| a2 /( a+1 )] 



e,- = e -6(j-i)(«+i)/«(i v iv^i^Vj). 



(3.15) 



(3.16) 



Proof. We only need to show ( 13.15ft and (13. 16ft hold when x — > oo along the sequence 
c n := n^ a+1 ^ a2 . By Proposition 13. 121 it follows that X ( a+1 >' a is regularly varying with the 
index a 2 /(a + 1). More precisely, as n — > oo, 



P(X { a+1)/a >x) 



^_ x -« 2 /(«+i)_ 



a 3 b 2 T(— a) 



(3.17) 
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By RemarkEHl we have E[\V k \ r ] < oo for any < r < a. Note that ||H fc || = x£+ 1)/a (lV 
V k ). By (13.171) and Breiman's Lemma, as n — > oo, 

PCIIHxH > x) ~ E[l V [T/ x r 2/(a+1) ] 3 " a " . x-" 3 /^; (3.18) 

ct J o z i (—a) 

see, e.g., Resnick (1987, p. 231). Then ||Hi|| is regularly varying with the index a 2 /(a + 1). 
By Lemma [37181 for any < r < a 3 /(a 2 + 1), we have E[||Hi - Hi|| r ] < oo. Now let 
H = (Hi, • • • , Hjfe) and H = (Hi, • • • , Hfc). By Markov's inequality and (13.181) we have, 
as x — > oo, 

P(||Hi - H x || > x) P(||H - H|| > x) 



V{\\K 1 \\>x) ~ P(||Hi||>x) 

< x I VE[||H 7 --H 7 -|rl ->0. 

~ P Hi >x) 4-f 1 J J J 

As in the proof of Lemma 3.12 in Jessen and Mikosch (2006), we obtain (13.151) from (13.181) . 
From the above relation we also have, as n — > oo, 

P^H G -HlHill > c n ) ~ P(c; 1 H G -HlHill > c n ). (3.19) 

Let H fc = (X e- f,(fc - 1) ) (Q+1)/a (l, V k ). For 5 > and K > 1 we have 

P(||H fc -H fc || >c„5|X >a n ) 

< P(||H fc - HfcH > cj, \V k \ < K\X > a n ) 
+ P(||H fc -H fc || >c n 6,\V k \ >K\X >a n ). 

Let J\ and J<i denote the two terms on the right-hand side. Then 

J x < V(K\{X e~ b ^)^/« -XtX 1)/a \ > c n 5\X > a^. 

By Theorem 13.151 and the continuous mapping theorem, we have J\ — > as n — > oo. Since 
Xq is independent of V k , we have 

J2 < P(|T4| > K\X > On) = P(\V k \ > K), 

which tends to zero as K — > oo. Then the regular variation property of X implies, for 

any £ > 0, 

lim P(||H fe - HfcH > x {a+1)/a C\X >x) = 0. (3.20) 

See, e.g., Resnick (1987, p.14) for a similar method. By (p£E]) we have P(K a ^ a+1 ^>X > 
a n) ~ h\TT x as n — y oo for some constant hi > 0. By (I3.20p and the multiplicative 
formula it is easy to see 

lim nP(||H fc - fifcH > c n 5, K a/(a+1) X > a n ) = 0. 

n— >oo 

It follows that 

limsupnP(||H fc - H fc || > c n 5, ||Hi|| > c n ) 

n— >oo 
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= limsupnP(||H fc - H fc || > c n 5,X { a+1)/a (l V V 1 ) > c r , 

n— >oo 

< limnP(X^ +1)/Q |^|l { | yi | >x} >c n ) 

= lim nEfl^r 2 /^)! ]P ( X ( a+1 )/ Q > c n ) 

n— >oo 

= c 1 E[\v 1 r 2 /^ 1 h UVll>K} ], 



where we have used Breiman's Lemma again for the second equality. The right hand side 
goes to zero as K — > oo. But, by (13 . 1 8[) there is a constant h 2 > so that P(||Hi|| > 
c n ) ~ h 2 n~ 1 as n — > oo. Then 

lim P(||H fe - fifcH > c^HlHxIl > c n ) = 0. 

n— >oo 

Let H = (Hi,-- - ,H fc ). We have 

k 

P(||H-H|| >c n «J|||Hi|| >c n )<X)P(||H i -H i || >c n 5|||H 1 || > c n ) -> 0. 

i=i 

Since Hi = Hi, by the above relation, we have as n — > oo, 

P^fi e -HlHill > c n ) ~ P^H G -HlHill > c n ) 

r-^nPdHHill >c„,c; 1 He-). (3.21) 

By Proposition 13.121 one can see, as n — > oo, 



o 



l+o 2 /(a+l) ' 



U 



Note that X is independent of for fc > 1. By the extended Breiman's Lemma, as 

72 — > OO, 

nP^H e •) — % E : «(E 1 , • - • , 3 fc ) e ■}] . 

where E,- = e - 6 (i-i)(°+i)/«(i ) ^) ; see Theorem 3.1 of Hult and Lindskog (2007). By (EDS]) 
and (13.211) . we have 

P^He -HlHill >c) 

^C 3 j( E[u{u > (1 V : u{B lt ■ ■ ■ , B») e •}] ^ 1+a2/(a+1) ■ 

Then (13 . 16[) follows by an application of the continuous mapping theorem. By Theo- 
rem [331 the sequence {H^} is jointly regular varying with index a. □ 



4 Point processes and partial sums 

In this section, we will first prove some limit theorems on the point processes associated 
with the stationary sequences and {H^} defined by (13.41) . From the limit theorems, we 
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derive the limits of suitably normalized partial sums for those sequences. The techniques 
have been developed extensively by Basrak and Segers (2009), Davis and Hsing (1995), 
Davis and Mikosch (1998) among others. For i.i.d. random variables, the idea goes back 
to Davis (1983) and LePage et al. (1981). Throughout this section, we assume 1 < a < 2. 
Let a n = n 1/a and c n = n (a!+1)/a2 = a { n +1)/a for n > 1. 

Lemma 4.1 Let r n = [n 5 ] for any < 5 < 1. Then for any x > we have 

lim limsupn > P(|£fc| > a nX, > a n x) = 0. (4.1) 

k=m 

Proof. Let E x denote the expectation of {X t : t > 0} given X = x. Take a constant 
r G (5, 1). For k >2, we can use Markov's inequality and Proposition 12.91 to see 



P(|£fc| > a n x, > a n x) 



a 



< 



[a n x)' 

(a n x) 
da 



-E 



L {|ei|>a n a;} 



e -b(k-s) X l/<* dZt 



k-1 



rot 
ret 



-E 



-{|ei|>a„x}Ex fc _ 



-E[ 1 {|ei|>a„x}(l + X T k ^ 



[CL n X 

£ (^ E [ 1 *-»-» (1 + A ' :e " < " 2>/<> 



where the last inequality follows from Proposition 12.101 In view of (II. 6p , we have 

r n 

n ^ p (kfcl > a nX, ki| > a n x) < Ji + J 2 + J 3 , 

where 



k=m 



J\ = ~, x — (r n — m + l)P(\ei\ > a n x), 



C 5 a ra n -rb(k-2)/a t?( yri \ 



k=rn 
oo 



k=rn 



By Theorem 13.161 we have P(|ei| > a n :r) = 0(n x ). It follows that Ji = 0(n <5 r ) as 
n oo. By Proposition 13.21 one can see, as n — > oo, 



E(\£i\ r l{\£i\>a n x}) ~ ^7— ;(a n x) r P(|£:i| > a n x) 

r ^P(\e 1 \> a n x). 



a — r 

1 

a — r 



Thus we have J2 = 0(n r ' a r ) as n — > 00. By Markov's inequality and Proposition 
E(X r l { | £l | >ana;} ) < 1 Q(1 _ r) E(Xg|ei| a(1 - r) ) 
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1 



— yEj^EUN^)]} 



< 



(a n x)«( 1 ~ r ) 

— -t: rE(Xn + Xo) 

(a n x) a ( 1 ~ r ) v u ; 



It follows that 



lim limsup J 3 < lim C 7 V e ~™Kk-i) = Q 

rn. — ^on „ , „ rn. — i(Yi < S 



k=m 



Then we have (14. ip . 



Theorem 4.2 Let G be defined by \3.12\l . Then we have, as n — > oo, 



□ 



(4.2) 



fc=l 



where t] is a point process on IRq tirai/i i/je Laplace functional E[e ''^j, / G C ( + ' 



o l K 6J # wen 



6/y 



exp 



oT(— a) 



E 



l-exp<^ - fly, 



y 



i + x 



G 



dy 

,,Q+1 ( ' 



(4.3) 



Proof. By Theorem 12.71 the process {X t } is strongly mixing with geometric rate. From 
(11.61) and (13. 4D we see is measurable with respect to a(Xk-i, X^). Then {1^} is also 
strongly mixing with geometric rate, and thus satisfies the mixing condition srf{a n ) de- 
scribed in Remark 13.91 with r n = [n s ] for any < 5 < 1. Since is a stationary 
sequence, we have 



nP( max \ek\ > a n x, \ei\ > a n x 

\m<\k\<r n 

— 71 [^d^l > anX ' > a " X ) ~^ ^*(\ £ -k\ > a n x , \ £ l\ > a n x ) 

k=m 

= n \P(\ £ k\ > a n x, > a n x) + P(|e 

k=m 
r-n+2 

< 2n P(l £ fc| > a n x i kl| > «n^)- 



k=m 



The right hand side tends to zero as n — > oo by Lemma 14.11 By Theorem 13.161 we have, 
as n — > oo, 



P( max |ejt| > a n x \eA > a n x) — > 0. 

\m<\k\<r n J 



By Theorem 13. 161 we have nP(||Ii|| > (cn) l l a ) — > 1 as n — > oo, where 

aa a (l -e~ ab ) 

c — . 

a 3 b 2 T(— a) 
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By Theorem 4.5 in Basrak and Segers (2009), we have (14.21) with the Laplace functional 
E[e~ v( -^} given by 



exp 



E 



E(G) 

This clearly coincides with (14.3 



l-exp<j - f[c 1/a v 



c l/a v 
1+Xr 



G 



d(—v~ 



□ 



Based on the above theorem, we now study the convergence of some partial sums 
associated with the sequence {1^} defined by (13. 4p . To do so, let us introduce some 
notation. For any B £ 38(M. + ) define 



n n 

U ltn (B) = 5> fc ls(M), U 2 , n (B) = E 1 + x. 1b ( 



£k 



1 + 



k=l k=l - ' Xk ^ 

Then we define U j>n {B) = U j>n (B) - E[U j>n (B)} for j = 1,2. 

Lemma 4.3 For any 5 > we have 

lim limsup P(a~ 1 |L r i jn (0, a n z] \ > 5) = 0. 

Proof. Since E(e fe ) = E(e fc |#fc_i) = 0, we have 

n 

^^1^(0, a n z\ = a" 1 ^[£fel{| efc |<a nZ } - E(e k l { \ £k \< anZ} )} 

k=l 
n 

= a ~n ^2[ £ kl{\e k \<a n z} - E(e fe l { | £fe \< anZ} | & h -\)\ 
k=l 

n 

~ ^2 \^( £ kh\£ k \>anz}\^k-l) ~ E(£fel{| £fc |> an2 }) 



(4.4) 



k=l 



Let Ji and J 2 denote the two terms on the right-hand side. By Theorem 13 . 1 61 one can see 
that e\ is regularly varying with index a/2. Then by Proposition 13.21 it follows that, as 
n — > oo, 

n 

Var(Ji) = a- 2 J2H[ £ khe k \<a nZ }-^kh\s k \<c^z}\^k~i)] 2 } 



k=l 



2-a 



< na n E(e 1 l { | ei |< an ^ } ) = na n E(£ 1 l {e 2< a 2 2 2 } ) ~ Cz 

which goes to zero as z — > 0. Now we discuss the asymptotics of J 2 . Observe that, for 
u > 72; + p, we have \Xi — ^x — p| > u if and only if X± > u + 7a; + p. It follows that 

u^Exd-Xi -ei|l{| B1 |> u }) = « a_1 Es(|7:r + p|l{|jri-^-p|>u}) 

= (72; + p)w a_1 P x (X 1 > u + 72 + p). 

Using Proposition 13.101 we see the right-hand side tends to zero uniformly in x £ [0, K] 
as u — > 00. By Proposition I3.11[ we have 

lim w a-1 Ej|£:i|l { | ei | >u} ] = lim u a ~ l E x [e^^-yx-p^u}] 
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lim u a 1 E x [X 1 li Xl>u+ y X+ p\} 

u— >oo 



(a-l)T(-a) 

and the convergence is uniform in x G [0, K\. It follows that, as n — > oo, we have almost 
surely 



1 1 {X k - 1 <K}~E[e k l{\e k \>a n z}\^k-l] 
n 

1 ^2 1 i x k-i<K} E Xk-i [ £ ll{|£i|>a„ 2 }] 



k=l 



k=l 



a a z l-a 



a 



{a _ 1)r( _ a)n E kx^<K}{qa + PA-l) + 0(1) 
!) r (_ a ) w zZ Vfc-i^}^ + Pa^fc-l) + °(1) 

( a _ i)r(-a) E [ 1{Xo -^ } ^ Q +^ x °)] + °W' 
where the last equality holds the ergodic theorem. Similarly, we have 
raa~ E[l{x <jr}£il{| ei |>on«}] 

= ( Q _ 1 ^ r ^_ Q ^ ) E [ 1 {Xo<A-}(ga+P a ^o)] +0(1). 

Then ( 14. 5ft and (14. 6ft cancel asymptotically as n — >• oo. Observe that 

P(1 {M} N > u) = P(a£ l {Xo>K} e- b ^Xy_ a dZ s > uj 

+ P(a^ 1 l {Xo>K} e- b ^xlL a d(-Z s ) 

By Remark \3M as w — > oo, 

P(l {Xo >jf>ki| > M ) ~ C ( K ) \P{Z X >u) + Y(-Z x > u)] = C(K)u~ a 

where 

C(K) = a a E\ [ l {Xo>K] e- ab{1 - s) X s ds\ < dE[l { x >x}(l + X )}. 
Then by Proposition 13.21 as n — > oo, 

n 

fl « lE { E [l{^/t-i>-ff} E (l £ fc|l{|£/tl>an2}l^fc-l) 

fc=l 

+ E ( 1 {^-i>A'}kfc|l{| £fc |>a„ Z })] } 



> It 



(4.5) 



(4.6) 
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2a n 1 E ( 1 { X k-i>K}\ £ k\ 1 {l { x k _ 1 >K}\e k \>a n z}) 



k=l 



a n 1 nE[l {Xo>Ky \e 1 \l { 



l {x >K}\ei\>a n z} 



= C 2 nP(l {Xo>jr} | £l | > a n z) = C 2 C(K)z~ a . 
The right hand side goes to zero as K — > oo. That gives the desired result. 

Lemma 4.4 For any 5 > we have 

lim limsup P (a~ 1 |f/ 2 ,n(0, a„2;]| > 5) = 0. 



□ 



Proof. It is simple to see that 



u a_1 E, 



x + 1 



-{|ei|>(a:+l>} 



[u{x + l)] a - 

(x + iy 



■E a .[eil{| ei | > ( a . + i) u }] ) 



where u(x + 1) > u and + a < 1. Thus as u — > oo, uniformly for x G [0, K], 



M a_1 E, 



x + 1 



L{|ei|>(x+l)u} 



0" a (g Q + p a x) 



(a - l)r(-a)(x+ 1)°' 



The remaining argument is similar to the proof of Lemma 14.31 



□ 



Theorem 4.5 Lei Z7i |n = Ui jn (0, oo) and [7 2in = £/2,n(0, oo). Taen we nave, as n — > oo, 

< 1 (^l,n,^2,n) A(C/ l5 C/ 2 ) On R 2 , 

where (£/i,£/ 2 ) iae a-stable random vector with characteristic function given by 
E[exp{i(AiC7 1 + A 2 Z7 2 )}] 

A 2 



exp < — E 



Ai + 



1 + X 

and p a = p a (l) and q a = q a (l) are defined by (COJ). 



(q a + p a X ) 



-iira/2 



(4.7) 



Proof. Fix z > and A = (Ai, A 2 ) G R 2 , and define the function on R 2 by /a,z(xi,x 2 ) 

AiXil{| :ri | > ^ + A 2 Z 2 l{|a!a|>*}- Tnen We have 

2 

Vn(fx,z) = a~ l ^2 ^jUj,n(a n z, oo). 

3=1 

It is easy to see that the mapping from M(R 2 ) into R defined by 
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fc=l 



is a.s. continuous with respect to the distribution of the limit point process r\ in Theo- 
rem EOJ By the continuous mapping theorem, as n — > oo, , we have r] n (j\ tZ ) r)(f\,z), 
and hence 



E 



exp|m n 1 ^A j L r J> (a n 2;, oo)| = E[exp{irj n (f x , z )}] ->■ E[exp(iv{f\,*))] ■ 
3=1 



where the right-hand side is given by 
1 



cxp 



oT(— a) 



E 



exp <^ i\iyl{ y>z y 



i^2yl{ y >z{i+x )} 



l + X 



-1)G 



dy 



By Proposition 13.21 and Theorem I3.16[ as n — > oo, 



a„ 1 E[i/i in (a n z,oo)] = na n 1 E(e 1 l { \ El]>anz} ) 

1 



anz . . 
P(|ei| > a n z) 



a (a — l)r(— a) 
1 



a — 1 
E{G)z 1 - a 

dy 



V(G)yl {y>z} 



aT(—a) J 

By Proposition 13.21 and Remarks 13.71 and 13. 8[ as n — > oo, 



a+l ' 



a n x E[U 2>n {a n z, oo)} = na n 1 E 

anz 
~ J 

a — 1 



If 



1 + X 
-E 



G 



a(a-l)r(-a) L(l + X 
1 



„l-a 



/■oo 


r G i 


/ E 


[(1 + X )«J 


/o 



1 ffy 

y l {y>*} „. a +l 



Consequently, as n —¥ oo, 



E 



z 

exp|m^ 1 ^A J f7 J> (a n z, oo)| = E[exp{i^ n (/ Aj2 )}] 



converges to 
exp 



aT(— a) 



i=i 



E 



cxp 



l\iyl{y>z} + j— y l{y> Z (l+X )} J 



-i\iyl {y>z} - - - — l{ y>z (i+x )} 
As z — > 0, the above quality tends to 



G 



J 2/ 



a+l 



exp 



1 



«r(— a 

= exp 



E 



i\iv+^- i\ 2 y 
e " 1+x o — 1 — ^Aiy 



aT(—a) J Q 



(e iz - 1 - iz)E 
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Ai + 



A 2 \ a 



1 + Xr 



dz 

yOt + l 



By Corollary 14.11 of Sato (1999) and ( 13. 3 j) one can see this coincides with (14. 7p . Since 
E(L r Jjn ) = E[Uj >n (0, oo)] = 0, by the above calculations and Lemmas 14.31 and I4.4[ as 

n — > oo, 



E 



exp |ia n i y^A i [/ J|W 



E 



exp jz'a^ 1 XjU jin (0, oo] 

3=1 



converges to (14.71) . That gives the desired result. 



□ 



Lemma 4.6 Let r n = [n s ] with < 5 < 1. T/ien we /icrae 

lim limsupnPf max ||Hfc|| > c n x,X > a n x) = 0. 

m— >oo n^-oo V — r n <k<— m / 

Proof. Since {(H/«,Xfc) : G Z} is a stationary sequence, by (jl.5p . it is easy to see 



nP( max ||Hfc|| > c n x, X > a n x 

V —r n <k<—m / 

= nP[ max ||HJ| > c n x,X m > a n x 

\m-r„<k<0 J 

< nP (e~ bm X + a J e- b{m ~ s) ds > a n x/2 ' 
+ nP[A n ,a J e- b ^XlL a dZ s >a n x/2j, 



where 



A n = ( 

By Theorem 13.191 it is easy to see that 



max IIHJI > c„x ) . 

r„<k<0 



lim P(A„) = lim VP(||H_ fc || > c, 

n— >co n— >oo 1 4 



lim (r„°+ 1)P(||H || > c n x) = 0. 

n— >oo 



Then Lemma 13.141 implies that 

limnP(A n ,a / e^-^ X]p 6Z a > a n x/2 ) = 0. 
From Proposition 13.121 it follows that 

lim nP(e- bm X + a [ e- b{m ~ s) ds > a n x/2] = Ce 
which goes to zero as m — > oo. Then we have the desired result. 



abm — a 



□ 
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Lemma 4.7 There exists 5 G (0, 1) so that for r n = [n s ] we have 

lim limsupnPI max ||Hfc|| > c n x,X > a n x) = 0. 

m— >oo n ->oo V m<k<r n / 

Proof. Recall that Xq is regularly varying with index a. It is easy to see that 

nPl max ||Hfc|| > c n x, X > a n x ) 

<n^ p (ll H fcll > c n x,X > a n x) 

k=m 

<n£P(||H fc -H fc || >c n x/2) 

k=m 

+ n^P(||H fe || > c n x/2,X > a n x). 

k=m 

Let Ji and J 2 denote the two terms on the right-hand side. We can choose r n = [n 5 ] for 
sufficiently small 5 G (0, 1), and thus 

lim sup ^limBup^EdlHi-Hiin 

n— >oo n—>oo X C n 

< lim sup — — E(||H 1 -H 1 |n = 0. 



n— >oo X r C r 



By Proposition 13.21 we have 



E[X l {Xo>an x}} ~ ^-yP(X > a n x) ~ C(a n x) 



l-a 



By Remark EH3 we have E[\l V V k \ a/{a+1] ] < 00. Note that (X ,X fe _i) is independent of 
Vfc for k > 2. Then for some constant 5 G (0,1/ a), 

2a/(a+l) 

J2 < /(a+1) £ 1 1 V ^r/^;X > a n x] 

k=m 
r i Tn 



-1) 



fc=m 

It follows that 



00 



lim sup J 2 < Cz l ~ a J2 e 



-b(k-l) 
lo—ri^j k=m 

which goes to zero as m — > 00. □ 
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Lemma 4.8 Let H fc = x£+ 1)/o (l, V k ). Then there exists 5 G (0, 1) so that for r n = [n s ] 
we have 

lim limsupnPI max ||Hfe|| > c n x, ||Hi|| > c n x) = 0. 

m->oo n-j.oo V m<\k\<r n J 

Proof. For any K > 1, we have 

nP max IIHJI > c n x, ||Hi|| > c n x ] 

Vm<|fc|<r n / 

< nP (X^ 1 '/ 4 * | VI 1 1 { | vi[>JO > °n x / 2 ) 

+ nP( max ||H fc || > c n x, KX { a+1)/a > c n x/2] 

\rn<k<r„ / 

+ nP( max ||H fc || > c n x, KX^ +1),a > c n x/2] . 

\ —r n <k<—m / 

Observe that X^ +l ^ a is regularly varying with index a 2 /{a + 1). By Remark 13.171 we 
have E[|Vi| 6 ] < oo for some b > a 2 /(a + 1). It follows from Breiman's Lemma that 

lim nV(X^ +1)/a \V 1 \l { \ Vl>K} > Cnx/2) 

n— >oo 

= lim nE^f'^h m>K} )P(X^ +1)/a > c n x/2) 

n— >oo 

= Cx -aV(«+l) E( | Fir 2 /(a+l) 1{| ^ |>x}) . 



The right-hand side goes to zero as if — > oo. Then the result follows by Lemmas 14.61 
and S3 □ 

Theorem 4.9 Let {Vj} be defined by \3.1$ - Then we have, as n — > oo, 

n 

^ : =E^H fe ^£ on M(Rg), (4.8) 
fc=i 

where £ is a point process on Mq with Laplace functional E[e~^], / G Cq~(]Rq) given by 

exp{ -f;/(y<° +1 >/'>e-'f- I K° + "/'»(l,V,))}]4| T }. (4.9) 

3=2 y ) 



E 



Proof. This proof is based on Theorem 4.5 of Basrak and Segers (2009) similarly as the 
proof of Theorem 14.21 It is easy to see that 

nP max ||H fe || > c n x, ||Hi|| > c n x) 
\m<\k\<r n J 

<nP(||Hi-Hi|| >c n x/2) 

+ nP( max ||H fc || > c n x, ||Hi|| > c n x/2 

V m<\k\<r n 
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By Lemma [3.181 we have E[||Hfc — Hfc|| r ] < oo for some r > a 2 /{a + 1). Then Markov's 
inequality implies that 

limsupnP(||Hi - Hi|| > c n x/2) < limsupn C - r E[||H 1 - H 1 || p ] = 0. 



By Lemma [4.81 we have 

lim limsupnPI max ||Hfc|| > c n x, ||Hi|| > c n x) = 0, 

where r n = [n s ] for some 5 G (0, 1). Let 



h 



aa 



-E[l V |T4| a2/(Q+1) ]. 



a% 2 Y(-a) 
By Theorem 13.191 we have, as n — > oo, 

nP{||Hi|| > (hnY a+1 ^ a2 } = nP{||Hi|| > /i (a+1)/a2 c n } -)• 1. 
Observe also that 

n n 

Uf) = E^ c n lH *) = Y,f( hia+1)/a2 ( hia+1)/a2c ^~ lH ^- 

k=l k=l 

Let ©j be defined as in Theorem 13.191 Then we can use Theorem 4.5 of Basrak and Segers 
(2009) to obtain g3J) with E[e _€ ^] given by 



cxp 



I f°° r r °° "i 

17^) 1 E[exp{-£/(^„e,)} 

3 — 2 

i-ex P { -/(/i (a+1)/a %ei)}Vi v |^r 2/(a+1) )]c/(-t;" a2/(a+1) ; 

— ^tt-TTT / E[exp{-y/(^ +1 )/ Q 4 )) 

E(l V \Vi\ a 2 /("+i)) J Q [ I ^ V J 



exp < — 



J'=2 



l-exp{ -/(u (Q+1)/Q 0i)})(l V |yir 2/(Q+1) )]d(-u- a ) 



exp 



i=2 



i-ex P {-/(J a+1 )/ Q e 1 )}Viv|y 1 | a2/(a+1) )]d(-0 

"5^r Efexpl-g/fe^vl^De,) 



exp 



3=2 



□ 



(l - exp { - f(y {a+1)/a (l V |T^ 1 |)© 1 )})]d(-y- Q! )|, 

which can be rewritten as (14. 9|) . 

From the above theorem, we can derive some limit theorem of partial sums associated 
with the sequence {H fc } defined by (13.41) . For B e =S^(R + ) define 

n n 

Si, n (B) = J2 X LiM x k-i), S 2 , n (B) = ^2x k - 1 e k l B (\X k . 1 e k \). (4.10) 



k=l 



k=l 
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Lemma 4.10 For any 5 > we have 

limlimsup P(c~ 2 \Si n (0, c n z) \ > 5) = 0. 



Proof. By Theorem 13.151 it is easy to see that Xq is regularly varying with index a/2 < 1 . 
Using Proposition 13.21 and Theorem 13.151 we have, as n — > oo, 

i n 2 

E[c£ 2 S liB ((W)] ^jE^iWiM] ~ ^P(Xo > c n z) ~ Cz 2 - Q 
The right-hand side tends to zero as z — > 0. Then we have the desired result. □ 

Lemma 4.11 Suppose that 1 < a < (1 + y/E)/2. T/ien /or an?/ 5 > we have 

limlimsup P(c~ 1 1 ^^(O, c n z) | > 5) — 0. 

Proof. By Theorem 13. 19} we see X £i is regularly varying with index a 2 / (a + 1). Under 
the condition 1 < a < (1 + v / 5)/2, we have a 2 /(a + 1) < 1. By Proposition 13.21 and 
Theorem 13.191 as n — )■ oo, 

1 n n 
E[c^ 1 |S , 2 , n (0,c n z)|] < — ^E[|X fc _i£: fc |l { | X;; _ iefc | <Cnz} ] = — E[|X £i|l { | Xo£1 | <CnZ} ] 



fc=l 
„2 



a + 1 — a 



2 " 



The right-hand side tends to zero as z — > 0. That gives the result; see also Davis and 
Hsing (1995, p.896). □ 



Theorem 4.12 Let V\ be defined by Let Si >n = Si >n (0, oo) and S 2t n = S2,n(0, oo) 



If 1 < a < (1 + y/E)/2, th en we have, as n — >■ oo, 

(a" 2 ^, c" 1 ^) — > (Si,S 2 ) on M 2 , 
where (S±, S 2 ) has characteristic function E[exp{iAiSi -MA2S2}] given 6y 



exp 



E 



exp 





/•oo 




7o E 


a ) 




ic 


" 2b A l2 / 2 


1 


_ e -2b 



(1 _ e -&(a+l))l/e 



y a+l 



(4.11) 



Proof. We first remark that the integral on the right-hand side of (14. lip is well-defined. 
In fact, by Remarks 13.71 and 13.171 we have, as x — > 00, 

P(Vi > x) ~ C lX ~ a + o(x" a ), P(Ui < -x) = o(aT Q ). 
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By Theorems 8.1.10 and 8.1.11 in Bingham et al. (1987), we have, as A — > 0, 

E[l - cos(AVi)] ~ C 2 \ a + o(A Q ), E[sin(AK)] ~ C 3 X a + o(A Q ). 

It follows that E(l — e lXVl ) ~ c\ a as A — > 0. Then the integral in (14. lip converges. Fix 
any A = (Ai, A2) € M 2 and z > 0, define the function on IR + x I by 

5^,2(^1,^2) = AiX^ a ^ a+1 ^l| x . 1>2 ( Q +i)/a} + \2X2^-{\x 2 \>z}- 

It is easy to check that 

£n(g\,x) '■= / g\,zd£ n = \ia~ 2 S hn {a n z, 00) + A 2 c~ J S 2t n{c n z, 00). 

JR+xR 

On the other hand, one can see the mapping from M(M.q) into R defined by 

00 „ 
N ■= E <W^) ^ N (9\ Z ) ■= / 9x,JN 



k=l 



is a.s. continuous with respect to distribution of the limit process £ in Theorem 14.91 By 



the continuous mapping theorem, as n — > 00, we have £, n (g\,z 

E[exp(i£ n (gx t ,))] -> E[exp(z£(sv))], 
where the right hand side is equal to 



£(&\,.z), and hence 



cxp 



Ejexp^A^^e" 26 ^ 1 )! 

3=2 



{ ye -6(i~i)> 2 } 



a 2 6 2 r(-a 

exp|zA 2 y v 2^ e • / i- L {| e - i '0-i)Cc«+i)/«^|> z } 

3=2 

E(l - exp{zA l2 / 2 + zA 22 /( a+1 )/ Q ^}^ ^ 



Then we can use dominated convergence theorem to see that, as z — > 0, 

ao" c 



cxp 



E 



a 2 b 2 T( 



— / E(l-exp{a 1 y 2 + a 2 ^ +1 )/^ 1 } 



exp (We" 26(i-1) + A 2 y( Q+1 )/ a e-^- 1 )( a+1 )/ Q \/,) 

3=2 



J y 



dy 
0+1 



Since the sequence {Vi, V 2 , ■ ■ ■ } is i.i.d., the above quantity is equal to (14. lip . Note that 
E(Vi) = 0. Then the theorem follows by Lemmas 14.101 and 14.111 □ 
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5 Asymptotics of the estimators 



In this section, we investigate the asymptotics of the estimators for the SCIR-model. The 
results are presented in a number of theorems. In fact, we shall first study the asymptotics 
of the estimators of the parameters (7, p) defined in (11 .4p . Their CLSEs can be obtained 
by minimizing the sum of squares in ( 11. 81) . They are given by 

- = Y2=i x k-i Y2=i x k-n Y!=i x k~i x k , 5 ^ 

(YJl=i X k-i) 2 - nY2=i x l-x 



and 



1 n n 

-[E**-7»E**-i]- (5 - 2) 



n 

fc=i fc=i 



By minimizing the weighted sum in (11.111) . we obtain the WCLSEs of the parameters: 

Ek=i(**-i + 1) Efc=i x^TTT ~ n 



and 

Pn 



n n 

-[E^-TnE^-r]- (5.4) 



n 

fc=i fc=i 



In view of Proposition 12.31 and the above expressions, in the discussions of the above 
estimators it suffices to consider a stationary realization {X t : t > 0} of the SCIR-model. 

Lemma 5.1 We have, as n — > 00, 

1 n 1 n 1 

-Y ^X, (5.5) 

k=i fc=i ft 1 



and 



1 n y 

i £^x7^pA + 7 (i-A) ( 5 - 6 ) 



where 



A = E 



1 + X r 



Proof. By Theorem l2.6l the process {^t} is exponentially ergodic and thus strong mixing, 
so the tail cr-algebra of the process is trivial; see, e.g., Durrett (1996, p. 351). Recall that 
E(-Xq) = a/b. In view of (II. 6p . we have 



El 



(rlkHKT^)+<fkH A ^ (1 - A) - 

Then the result follows by Birkhoff's ergodic theorem; see, e.g., Durrett (1996, p. 341). □ 



35 



Theorem 5.2 The estimators (p n ,ln) are strongly consistent and, asn — >■ oo ; n( a 1 ^"(7 n " 
7, p n — p) converges in distribution to 

i 71 ~ 1 (f/i,f/ 2 )( \ ^"i 1 ) =F' 1 (XU 1 -U 2 ,(\-l)U 1 + ab- l U 2 ), 

where F — (1 + a6 -1 )A — 1 and (L^, C/ 2 ) an a-stable random vector with characteristic 
function given by fiJTfy . 

Proof. We first remark that (14.71) defines a Gaussian random vector (Ui, U2) when a = 2. 
In view of (15.31) and (15.41) . the results of Lemma [5.11 imply that 

a .s. ab~ 1 \ — pX — 7(1 - A) _ ab~ 1 'j\ - 7(1 - A) 
7n ~^ {l + ab-^X- 1 ~ (l + ao-^A- 1 ~ 7 



and 



Pn ^> |(1 -7) = P- 



Those give the strong consistency of p n and 7„. Write 

_ ELi(7^-i + 7)ELi^diTi-"ELi : feff 

Z^k=l l^k-l l^k=l + 1 lb l^k=l X k _ L +l 

ELi(**-i + 1) ELi xTTTT - 

ELi(7^-i + p) ELi x^n - ; ELi fe^ 

ELi(^ + i)ELi3d+i- n2 
Efc=i(^ ~ £ k) Efe=i x fe _!+i - n Efc=i Xk-1+1 

ELi(^i + i)ELix^tt-^ 2 ' 

where the last equality follows from ( II. 6p . Combining this with (I5.3P we have 



Pn - P 



n 
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(5.7) 



V n I V" £1.-71 V 1 — ^ 

Efc=i(**-i + !) E fc =i x^TTT - 72 
By ( II. 6p and (15. 4p it is easy to see that 

- _ ELl( X fc - l X k-l) - (in ~ 7) ELi X k-1 
Pn — 

n 

_ ELl(P + £ k) - (in - 7) ELi X k~l 

n 

Using (15. 8p we have 

ELi £ k - (in ~ 7) ELi X k-1 



En e k sr^n y _ sr^n sr^n X k _ 1 

k=l Z^fc=l ^k-l Z^fc=l fc fc Z^fc=l Xfc.j + l 



=1 

We can rewrite (15. 8ft and ( 15. 9ft into the matrix form as 

(in -l,Pn-p) = f„ _1 UA, 



rr 



where 



(a-l)/a 



1 n n -. 

fc=i jt=i s 1 



(5.9) 



B, 



^ ( 2^ = 1 Z^fc = l X fe _! + 1 

-n ELi 



and 



n n 



v fc=: 

From (15. 5p it follows that F„ F and 



k=l k=l K 1 



B„ ^> B 



A A - 1 
-1 ab~ l 



In the case 1 < a < 2, we have U n — (C/i, £/ 2 ) by Theorem 14.51 In the case a = 2, we 
have 



E[e?] = a 2 E 



aa2 rl e-^ds<^. 



Then, for any 5 > 
1 



n 



fc=i 



E 



Ai + 



a 2 



en. 



tends to zero as n — > 00. Furthermore, by the ergodicity theorem, as n — > 00, 



\ 2 _ 2 



fc=i 



(T 



+ 1 

Ai + 



A 2 ^ 



(pa + <?2^o) 



1 + X , 

Then a martingale convergence theorem implies that, as n — > 00, 

A 2 

fc=i 
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(5.10) 



1 " 



2 A' 



converges in distribution to a Gaussian random variable with mean zero and variance 
given by the right-hand side of (I5.10P ; see, e.g., Durrett (1996, p. 417). It follows that 



E 



i n A 



converges to the right-hand side of (14.71) with a = 2. Since (Ai, A2) G K 2 can be arbitrary 
in the above, we also conclude U n (Ui, 1/2). That proves the desired convergence. □ 

Theorem 5.3 The estimators (b n , a n ) are strongly consistent and asn — >■ oo ; n( a ~ x >l a (b n — 
b, a n — a) converges in distribution to 

F- 1 {e\U 2 - AC/i), (1 - e- 6 )- 1 ^ + b{\ - l)]U x + ab^e\U 2 - Atfi)), 

where F — (1 + a& -1 )A — 1 and (U±, U2) is an a-stable random vector with characteristic 
function given by ([^.7] ). 

Proof. The strong consistency of b n and a n follows from that of p n and 7„. By the 
relations in (11.4)) . we have, as n — > 00, 



(in - 7) = e" fe " - e~ b = -(b n - b)e~ b + o(b n - b) (5.11) 



and 



a n - a 



Pnbn pb p n 6 n (l - e b ) - pb(l - e 



\- e ~K l-e~ b (1 - e - 6 ™)(l - e" 

b n (pn ~ p) p(b n - b) , pb(e~ b ™ - e~ b ) 



l- e -bn l- e - fc " (1 - e- fe ")(l - e~ b ) 
KPn - P) ae b (%-i) a(7n-7) , ,r i \ ftl 1Q x 

= T3T ^ r 5 + T^r + o(6 " " 6) - (5 - 12) 

Then the desired convergence follows from Theorem I5.2I □ 

Theorem 5.4 The estimators (p n ,7n) ore weakly consistent. Moreover, if 1 < a < 
(1 + \fb)/2, then, as n — > oo, 

n {Q " 1)/a2 (7n-7,Pn-p) -^5^2(1, -aft" 1 ), (5.13) 
where (Si, £2) /ias characteristic function given by 

Proof. By (jl.6p and (15. ip we have 

. _ = ELi X k-i ELi ( g fc + P) ~ w ELi x fc-i( g fc + P) 

ln 7 „ ( ELi ^fc-i) 2 - n ELi x l-i 
_ ELi -^"fc-i ELi g fc ~ n ELi ^-i g fc 

(£Li**-i) 2 -» £Z=i*2-i 
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Then using (15. 2p we get 

* _ _ ELi X fc - 7n ELi _ . 

Pn P — P 

n 

= Tl=l( X k ~ T^fc-l) - (7n ~ 7) ELi X k-l 

n 

= ELi( g fc + p) - (7» - 7) ELi x fc-i 
= ELi g fc - (7n - 7) ELi x k-i 

n 

fc=i A fc-i 2^fc=i A fc-i g fc - 2^fc=i A fc-i Z^fc=i g fc 

(ELi^-if-^ELA 2 -! 

By Theorem 14.121 it is easy to see that 

ELi x k-i£k p 



En y2 
k=l ^k-l 



0. 



Then we have p n — p and j n — 7 0, giving the weak consistency of (p n ,j n ) 
From the above relations it follows that 



1 y->n y r 1 v^n y 

(a-l)/a 2 ^ _ ^ _ n i+(^+i)/^ Z^fc=l ^fc-l 2^k=l £ k g^+g/g 2-^fc=l ^-fc-l£fc 

„i+2/c« ( Ea,-i ^fc-i) — -hs Efc=i ^i-i 



n^" 1 ^ ( 7n - 7 ) = 

wn ,; 1 / y^n y ^2 ! ^ n 

n i+2/<* I, Z^fc=l ^fe-l, 

and 



1 I v^»t y NT^n y _ _ v^n y 2 v^n _ \ 

(a-lj/o 3 /* _ n ni+("+i>/° 2 I ^ fc=1 fc ~ 1 ^ fc=1 Afc ~ l£fc 2^fe=l A fc-1 2^fc=l £ fcj 

1 f Y^" y ~\ 2 _ 1 y2 
n l + 2/a ^ Z^fc=l n 2/a l^k=l^k-l 

Take any constant < 5 < [1 A (a — l) 2 ]/a 2 . We can rewrite the above relations into the 
matrix form 



(in -J,Pn-p) = T n 1 S n A T 



where 



\ n 2 \ n 



k=l k=l 



1 V 1 V^ra y-2 

/ Q 2_ a 2^fc = l ^fc-l „1 + 1/ Q 2_ 4 Z^fc = l ^fc-l 



A — / „l + l/^-5 ^fc = l"-K--l „1 + 1/ 

■**-TJ. I - 1 



fe-i 



and 



/ 1 - 1 - 

Sn = ( n l/a+8 ^2 £fc ' n (a+l)/a 2 

^ k=l k=l 



X k-l £ k ■ 



If a 2 < a + 1, by (j5.5p and Theorem |4.12| we have T n — > —Si and 

A AA-f ° ° 

n y — 1 afe" 1 

By Theorems S3] and E2 we have S n -A (0, S 2 ). Then fl5TT3|) holds. □ 
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Theorem 5.5 The estimators (b n , a n ) are weakly consistent. Moreover, if 1 < a < 
(1 + VE)/2, then, as n — > oo, 



n 



(«-D/« 2 (p n _ b ^ n -a)-^ -e\l, ab~ 1 )S^ 1 S 2 , 



where (Si, S 2 ) the characteristic function given by ^4- 11]) . 

Proof. The weak consistency of b n and a n follows from that of p n and j n . The relations 
(15. lip and ( 15 . 1 2 j) still hold when the "checks" are replaced by "hats" . Then the desired 
result follows from Theorem 15.41 □ 



Theorem 5.6 The estimator a n for a defined by fa!.14\ ) is weakly consistent. 

Proof. For p e (0, a) and 8 e (0, min{l — I/a, l/a 2 }), we need to show that, as n — » oo, 



ii 



J; I X t k - ^f fc -i 

l-p/a | yl/a , 8 

k=i A t k -i -r n 



A o-"E[|Z 1 | p ]. 



(5.14) 



By equation (jl.2|) . we have 



l.h.s. of (EH 



— Y 

k=l 

n 

1 V 

„ l-p/a 



X t k - X t k -i 



l/a 



k=i ^t fc _i + n 



k=l 



Ik 



aX 



l/a 



4-" + n ~ & 



*** crl 



l/a 



n l-p/a jL^i 



k=\ 



l/a 



I /" 

+ — ; — 7- > / crdZ, 



tit 



p 


pt k 






/ <7<iZ s 











It is easy to see that 



1 n 

n l-p/a Z—r- 

fc=l 



l/a. 



1 " I —b /'* X s ds + o(t fc - t k - 

< 7- > 

r l—p/a Z-~/ 



X 1/a +n~$ 



< n p/a (b sup X, 
^ se[o,i] 



p 

a I n 



which goes to zero a.s. as n — > oo. Observe that 



n 

1 V" 
n l-p/a Z-^t 

k=\ 



t k 



aX 



l/a 



p 


rt k 







l/a 



^1/a 
L *fe-i 



10 



k=l "^fc-i 

By equation (jl.2p and Lemma 12.81 



E 



sup {Xg-Xt^l 
se[t k _ 1 ,t k ] 













< E 




+ - + E 


sup 








n 


L se[t fc _i,t fe ] 





ex xj_ dZ u 



<C 2 



sup X s 
se[o,i] 
1 1 



CicxE 



X s ds 



l/a-i 



n n 



l/a 



Then by Lemma 12.8 



n l-p5-p/a 
< 



E [E| / <xZ a - X&dZ, 



k=l 



1—pS—p/a 



n 



5> 



k=l 



tk-i 



X s -X t \ds 



fc=l ol k-l 



p/ai 



p/a 



pS—p/a 2 



<C 3 n 

The right-hand side tends to zero as n — > oo. It follows that J2 — as n 00. Using 
the stationarity of {Xj} and the self-similarity of {Z t } we have 



-rV E E 

n l-p/a I Z-^t 



n , rt 



n 



vl/a . 
1 *t*-i + « 



-dZ. s 



fe=i 



E 



n 



LVX 1/a + n^ 



which goes to zero by dominated convergence theorem. Observe that 



1/3 : -,1-p/a E 



1/a 



n x-p/, 

fc=l 

,1— p/a 



7T 



dZ s 

n- & 



p 


rt k 


V 










Jtk-i 





1 X tk _ x + " 



Then J3 — s> as n — >■ 00. On the other hand, it is not hard to see that 

n 1 n 

^ *£o*\z tk - z tk j p = - E^* - z *-tf ° p n\zi\ p \ 

k=l k=l 

Thus we have (15.141) . 



n 1 



□ 
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